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ABSTRACT 



We study the axion-type curvaton model, with emphasis on the large field 
regime where analytic results are very difficult to obtain. We evaluate the 
tensor-scalar ratio r using WMAP normalization and the non-linearity pa- 
rameters /tvl and g^L by solving the equations numerically using the 5N 
formalism. We compare them with results for the curvaton with quadratic 
potential. We find that r is much smaller for the axion-type case compared 
to the result from the quadratic potential, with the difference increasingly 
more pronounced at larger field values. gNL is found to be positive, unlike the 
quadratic case where it is negative, and the amplitude of gNi is much larger. 
Moreover, there is a nearly linear scaling between g^L and /tvl, with small 
deviation from linearity at large field values. The slope between gjyj^ and f^L 
depends on the parameters characterizing the axion-type curvaton model. We 
further consider a mixed scenario where both the infiaton and the curvaton 
contribute to the primordial power spectrum and the non-Gaussianity pa- 
rameters are found to be much larger than those in the case with quadratic 
potential. 



1 Introduction 



Despite the simplicity of the inflationary paradigm [1], the mechanism by which curvature 
perturbations are generated is not yet fully established. An alternative to the standard 
scenario where the inflaton fleld is responsible for the accelerated expansion of the universe 
as well as the generation of curvature perturbations, is the curvaton mechanism [2] . This 
scenario has been gaining increasing popularity recently. According to this scenario the 
flnal perturbations are produced from initial isocurvature perturbation associated with 
quantum fluctuations of a light scalar fleld, the so-called curvaton, whose energy density 
is negligible during inflation. These curvaton isocurvature perturbations are transformed 
into adiabatic ones when the curvaton decays into radiation long after the end of inflation. 
It is also conceivable that the flnal curvature perturbations arise from a mixture of inflaton 
and curvaton perturbations. The curvaton scenario may naturally lead to much higher 
levels of local-type non-Gaussianity than standard single-fleld inflation. 

The local-type non-Gaussianity can be characterized by some non-linearity param- 
eters, namely, J'nl, Qnl and so on. They appear in the non-linear expansion of the 
curvature perturbation as: 

C(X) = C,(X) + \fNLeM) + ^9nlCI{?^) + (1.1) 

where C, is the curvature perturbation in the uniform density slice, C,g denotes the Gaussian 
part of The primordial power spectrum, bispectrum and trispectrum are deflned by 

(C(ki)C(k2)) = {27:fVc{h)5\\^^ + k2), (1.2) 

(C(ki)C(k2)C(k3)) = {^TxfB^ih, k,)6'ik, + k2 + ks), (1.3) 

(C(ki)C(k2)C(k3)C(k4)) = i2nfT^ik,, k^, k^, k,)5-'{k, + k^ + ks + k4). (1.4) 

The bispectrum and trispectrum are respectively related to the power spectrum by 

Bciki, k2, k-s) = pNL[Vc{ki)Vcik2) + 2 perms], (1.5) 

Tc(/Ci,/C2,fc3,/C4) = rjvL[7^c(A;i3)Pc(^3)^c(^4) + 11 perms] 
54 

+ —gNL[Vc{k2)Vc{ks)V^{k^) + 3 perms]. (1.6) 
The parameter tnl depends on f^L, as 

if the total curvature perturbation is generated by only one fleld. The parameters Jml, 
TisiL and qnl parameterize the non-Gaussianty from the irreducible three-point and four- 
point correlation functions respectively. See [3] for more discussions on the bispectrum 
and trispectrum. 
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In the simplest version of inflation model Jnl ~ 0{ns — 1) [4]. It is constrained 
by WMAP5 (n, = 0.96l[J:[!J^) [5] to be much less than unity. Even though a Gaussian 
distribution of the primordial curvature perturbation is still consistent with WMAP5 
(—9 < /nl < 111) at 95% C.L.), the limit on fiy^ on the negative side has been greatly 
reduced from WMAP3 to WMAP5, while still allowing for large positive value. The latest 
limit on this parameter is Jnl = 38 ± 21 at 68% CL for the local shape non-Gaussianity 
in [6]. This has prompted many authors to consider different possible mechanisms for 
generating a large positive non-Gaussianity [7-13]. 

The curvaton is a light scalar field whose mass is small compared to the Hubble pa- 
rameter H^, during inflation. It is sub-dominant during inflation and its fluctuations are 
initially of isocurvature type. After the end of inflation it is supposed to completely 
decay into thermalized radiation before primordial nucleosynthesis and the isocurvature 
perturbations generated by the curvaton are converted into final adiabatic perturbations. 
In order to prevent large quantum corrections to the curvaton mass and keep it small 
enough during inflation some symmetries are called for. It seems natural to invoke su- 
persymmetry. However, supersymmetry must be broken above the inflation scale and 
the mass square of each scalar field generically receives a correction of order H^. So 
supersymmetry cannot be invoked to protect the smallness of curvaton mass. 

A promising curvaton candidate is the pseudo-Nambu-Goldstone boson - axion [14], 
whose potential for the canonical field a takes the form 

V{a) = m'f {1- cos j), (1.8) 

where / is called the axion decay constant. Actually the axion field is generic in various 
string models [15]. The smallness of the axion mass is protected by the shift symmetry 
a ^ a + 6. If cr <^ /, the curvaton potential is approximately quadratic, and the second 
order non-linearity parameter gisn is expected to be small. However, (^jvl is large if the 
non-quadratic term in the curvaton potential plays a significant role [9-11]. 

In the axion-type curvaton model, it is very hard to get analytic results if a ~ /. This 
regime of large field value was studied numerically in [12] where the authors calculated 
Jnl- They concluded that it is quite similar to that of the quadratic potential if the 
curvaton accounts for the observed curvature perturbations. In this paper the curvature 
perturbation in the axion-type curvaton model will be studied in more detail, and we 
will see that the WMAP normalization and the second order non-linearity parameter 
qnl are quite different from those of the quadratic one, in the large field regime. The 
curvature perturbation is calculated numerically using the SN formalism [16] and the 
results are compared with the quadratic model. We find that, as expected, r agrees 
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with the quadratic case at small field values. However, it is much smaller for larger field 
values and the difference increases rapidly as the field value increases. Qnl is positive and 
its amplitude is much larger than the quadratic case. We find that qnl scales linearly 
with Jnl, with small deviation from linearity at large field values. The linear behavior is 
similar to the quadratic potential but the slope is quite different. We also analyze a mixed 
scenario where the primordial power spectrum has contributions from both the inflaton 
and the curvaton. We find that giyL and f^i are greatly suppressed with respect to the 
case where only the curvaton contributes. However, the suppression is much less for axion 
compared to the quadratic potential. 

This paper is organized as follows: in Sec. 2, some analytic results for the curvaton 
model will be briefiy reviewed. In Sec. 3, we numerically calculate the curvature pertur- 
bation for the axion-type curvaton model based on the 6N formalism [16]. We calculate 
Jnl and qnl and demonstrate the results in comparison to the quadratic potential and 
the potential with quartic correction. Next, we discuss the mixed scenario where the 
total primordial power spectrum has contributions from both the infiaton as well as the 
curvaton. Our results are summarized in Sec. 4. In an appendix we present a general 
discussion for the curvaton model with a small correction to the quadratic potential. 



2 Analytic results for the curvaton model 

We begin with a brief review of the curvaton model with quadratic potential. In the limit 
(7 ^ / the axion-type curvaton potential is reduced to this case. The amplitude of the 
primordial power spectrum generated by the curvaton with quadratic potential 

V{a) = ^-m'a' (2.1) 



is given by 



1 H"^ 

-^C- ^ TT^fo^^ (2-2) 
'dTx^ at 



where H,, and a,, are the Hubble parameter and the vacuum expectation value of curvaton 
field evaluated at Hubble exit during infiation, respectively, and 

Id = -. 7=r—^ ^a,D = — (2.3) 

4-i2<x,D Ptot 

evaluated at the time of curvaton decay. WMAP normalization [5] requires 

P^,obs = 2.457l°;°^2 ^ ^Q-9_ ^2.4) 
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The amplitude of the tensor perturbation only depends on the scale of inflation, as, 

= (2.5) 

The tensor-scalar ratio is defined as r = Pt/Pq and then the WMAP normalization 
(P(^ = P(^^obs) implies that the Hubble scale during inflation is related to r by 

^ 1.1 X 10-V5Mp. (2.6) 

Usually the dimensionless parameter r is used to characterize the inflation scale H^, . If the 
total amplitude of the primordial power spectrum is generated by the curvaton, namely 
P^a = P^obs, then r is given by 

r = ^^ (27) 

which encodes the WMAP normalization. Since the potential takes the quadratic form, 
the curvaton field evolves linearly after inflation and the non-Gaussianity parameters /tvl 
and qnl are respectively given by [17], 

. _ 5 5 5/i5 

- 4^ " 3 " "T' ^^-^^ 

9NL = ^ [-j^ + 2 + + 3/^ j . (2.9) 
For large positive /atl, fo <^ 1, which implies 

g^L ^ -y/iVL. (2.10) 

Since qnl is the coefficient of Cg, it is very difficult to measure qnl in the near future if 
it is the same order of fNL- 

After inflation the universe is dominated by radiation and hence the Hubble parameter 
goes like a~^. Roughly, the curvaton does not evolve until H = m. Once the Hubble 
parameter drops below the curvaton mass, the curvaton field rolls down its potential and 
starts to oscillate around its minimum. The energy density of an oscillating curvaton field 
with quadratic potential goes like dust-like matter (~ a~^) and grows with respective 
to that of radiation. When the Hubble parameter is roughly the same as the curvaton 
decay rate To-, the curvaton starts to decay into radiation. For simplicity, the curvaton 
is assumed to suddenly decay into radiation at the time of if = Fg-. Let us assume the 
scale factor equals to 1 at the time of H = m, and the scale factor is denoted by a when 
H = Fg-. If a^^, <^ Mp, the radiation and curvaton energy densities are respectively 3Mpm^ 
and ^m'^al at the time of H = m, and thus. 



3M^m^a-'^ + ^mV.^a"^ = 3Mp^F^. (2.11) 



The curvaton energy density parameter at the time of its decay can be written by 



^^,n = ^a-\ (2.12) 



6M2r2 



If ^ < then, a ~ ( gT/jjr^ ) ^"^^ and ^I^^d - 1, which gives, /nl - -5/6. If, on the 
other hand, ^ ^ then, a = a/ m/V^, and 



This gives 



which is much smaller than 1, and hence we get /nl— 4^ ^ 1. In this paper we focus 
on the latter case with large positive Jnl- Combining Eqs.f l2.7p and (12.141) we get 



r 

m al 



r = 1152— (2.15) 



fNL and qnl are respectively given by 



/^^^10^/i^^, (2.16) 
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100 /r^ 3^ 

3 V ?n cr. 



g^^C^- — ^!^^. (2.17) 



We see that WMAP normalization implies 

^^^lO-^f-^-V. (2.18) 
m \JnlJ 

Then, if r ~ 10"^ and /jvl ~ 10^, Vjm ~ IQ-^^^ 

In order to compare with the results of the axion-type curvaton model, it is better to 
re-scale the curvaton field to be 6* = a//, and we get 



8 M2 V r„ 

In this paper, we choose f /Mp = 5 x 10~^ and T„/m = 10~^ which are the same as 
those in [12]. This gives, fo = 3.125 x 10"^ 9l, and r = 4.6 x 10^ 9~'^. For 9, < 1, 
we get r > 10'^, which is not compatible with the bound from WMAP5 of r < 0.20. 
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On the other hand, if r is required to satisfy the bound from WMAP5, the curvature 
perturbation generated by curvaton will be much smaller than Pc^^obs for 0* ^ 1. We will 
more carefully discuss this case in Sec. 3.2. If we tune Tcr/m to be very small, for example 
10-^\ r < 0.2 and P^^ = P^^ 

obs can be achieved. But numerical calculation becomes much 
more difficult as we decrease T„/m and is not much more instructive. So, in this paper, 
we keep f/Mp = 5 x 10"^ and Tu/m = 10"^ fixed and illustrate the physics. 
In [10], the cosine-type potential of curvaton is expanded to the second order. 



\/^-mT^^(^l--0^ + ...J, (2.21) 

and the second term in the bracket is taken as a correction for small 6 values. Taking 
into account the non-linear evolution of curvaton after infiation but prior to oscillation, 
the non-linearity parameters are given by 

fNL = ^{l + h)-l-^^, (2.22) 
4/d 3 6 

25 / 9 „, , 9 , , 1 



where 



9NL = ^ [jf^ih + 3/^2) - -^(1 + h2) + -(1 - 9/12) + lOfD + 3/^ ) . (2.23) 



1.1355(0.951 + 0.1895) 
(0.951 + 0.5685)2 ' 
1.1355(0.951 + 0.1895)2 



(2.24) 
(2.25) 



(0.951 + 0.5685)3 ' 
5 = ^- (2.26) 

If 6** ^ 1, the above results are the same as those of quadratic potential. Since /i2 is 
positive, fNL is enhanced due to the correction. Another important feature is that Qnl is 
expected to be large and positive for large 6^. Even though this expansion is not reliable 
for large values of 6 (see Fig. 1), (note, in particular, that the potential 6*^/2 — 6'^/24 
approaches its maximum value at around 6 = ±0.787r), we will naively use the above 
formula to calculate Jml and Qnl, and compare them to the numerical results in Sec. 3.1. 

Before we end this section, we would like to mention the vacuum expectation value 
(VEV) of the curvaton during infiation. At the classical level, the VEV of the curvaton 
is taken as a free parameter. However, in a long-lived inflationary universe, the long 
wavelength modes of the quantum fluctuation of a light scalar fleld may play a crucial 
role in its behavior, because its Compton wavelength is large compared to the Hubble 
size during inflation. In [18], the vacuum expectation value of the square of such a light 
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Figure 1: The red (solid), blue (dotted) and green (dashed) lines correspond to the axion 
potential, quadratic potential and 6*^/2 — 6^/24 respectively. 



where is the Hubble parameter during inflation. Typically the vacuum expectation 
value of a can be estimated as 



For a curvaton field with quadratic potential, its typical vacuum expectation value is 
given by Hl/m and its energy density is roughly H^. Larger the Hubble parameter 
during infiation, larger the energy density of curvaton. However, the energy density of 
axion-type curvaton is bounded by m? p from above. If mf ^ H^, a* can be estimated 
as Hl/m. On the other hand, a* ~ / if mf < H^. Since it is difficult to get analytic 
results for the axion-type curvaton model if a* ~ /, we will use numerical method to 
calculate the curvature perturbation in the next section. 




(2.27) 




(2.28) 
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3 Numerical results of the curvature perturbation 



The curvature perturbation on sufficiently large scales on the uniform density slicing can 
be computed by the 6N formalism [16]. Starting from any initial flat slice at time tmi, on 
the uniform density slicing, the curvature perturbation is given by 

C{t,^) = 6N = N{t,^)-No{t), (3.1) 

where N{t, x) = ln(a(t, x)/a(tj„j)), and No{t) = ln(a(t)/a(ti„j)) is the unperturbed amount 
of expansion. In the following subsections we will study two different curvaton scenarios, 
flrst, the curvature perturbation is generated by the curvaton alone, and secondly, the 
curvature perturbation is a mixture of inflaton and curvaton perturbations. 



3.1 Primordial power spectrum mainly generated by curvaton 

In this subsection we focus on the curvaton model in which the curvature perturbation is 
generated by a curvaton fleld. The curvature perturbation can be expanded as 

C = NJcr + ^N^^Ja^ + ^N^^^^ida)^ + ... (3.2) 

where N^^ = dN/da, N^^^ = d'^N/da'^ and so on. Since 5a = the amplitude of the 
curvature perturbation generated by curvaton takes the form 

Pc,-iVJ(§)^ (3.3) 
and the non-Gaussianity parameters are given by 

6 AT/' 54 



J^L - ^TTT, 9nl - TT^rr^- 



After inflation, the vacuum energy which governs the dynamics of the inflaton is 
converted to radiation and the equations of motion of the universe are 

H' = ^{pr + p.), (3.5) 

Pr + AHpr = 0, (3.6) 

P. = la' + Via), (3.7) 

a + 3Ha + ^ = 0, (3.8) 
da 

where pr and p^ are the energy densities of radiation and curvaton. It is convenient to 
deflne new dimensionless quantities, namely N{t) = In a{t), x = mt and 9 = a/ f. Then 



8 



the above equations of motion are simphfied to be 



dN 
dx 



(3.9) 



where 



and 



d'^d ^ ^dN de ^ dV{9) _ ^ 
dx'^ dx dx dO ' 



V{a) = m'rV{e), (3.11) 



Pr,ini r>\ 

"=3iW>^' ^'-''^ 
and pr,mi is the radiation energy density at t = tmj. For the axion-type curvaton model, 
V^(^) = 1 — cos^. The scale factor can be rescaled to satisfy a{tini) = 1 and then N{tini) = 
0. Roughly speaking, the number of e-folds from the end of inflation to the beginning of 
the curvaton oscillations is almost unperturbed and the time when H = = m can 
be taken as the initial time for simplicity. Because the energy density of radiation is 
still much larger than that of curvaton when H = m, Pr,ini = ^M'^m? and then a = 1. 
Here we adopt the curvaton sudden-decay approximation which says that the curvaton is 
proposed to suddenly decay into radiation at the time oi H = To-. This condition is used 
to stop the evolution of curvaton in the numerical calculation. We must mention that all 
our numerical results presented in the following are for / = 5 x 10~^Mp and T„/m = 10~^, 
unless otherwise specified. 

First of all, we want to compare the analytic results for the curvaton model with 
quadratic potential to the numerical results. In this case, the non-Gaussianity parameters 
only depend on the parameter We numerically calculate exact f£, which is given in 
Eq. (12. 3p and show it along with the analytic approximate expression for fo given in 
Eq. (I2.19P in Fig. [2l We see that the exact numerical result (blue dashed line) is larger 
than the analytic approximation (black dotted line) by about factor of two. The reason 
is that the energy density of the curvaton, as it oscillates about its minimum potential, 
does not exactly decrease as fast as a~^, due to the friction term 3Ha, even though this 
term is small compared to the driving term m^cr. So the curvaton energy density at the 
time of its decay is slightly larger than the analytic approximation. 

Using the dimensionless quantities used in the numerical calculation, the amplitude of 
the power spectrum of the curvature perturbation generated by the curvaton is denoted 
by 

n.=A^fJ^, (3.13) 
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Figure 2: The blue (dashed) hue shows the exact value of fo computed numerically, 
whereas, the black (dotted) line shows to fo computed from the analytic approximation 
in Eg. (12. 191) . for the quadratic potential. 



where 



N' 



tot 



dOir. 



(3.14) 



and Ntot denotes the total number of e-folds from tini to the time of H = To-. The 
tensor-scalar ratio r (or the inflation scale i/*) is determined by WMAP normalization as 

8 p 



(3.15) 



The numerical result for r is shown in Fig. [3] where we have plotted for axion-type (red 
solid lines) and quadratic case (blue dotted lines). The upper set of red and blue lines 
correspond to Tcr/m = 10^^, while the lower set corresponds to Tcr/m = 10~^. From 
this figure we see agreement of the value of r between the two models for 6^, <^ 1. As 6^ 
increases we find that r for the axion-type potential increasingly becomes smaller than 
the quadratic potential, r seems to scale linearly with T„/m for the axion-type model, as 
is the case for the quadratic potential, as seen in Eq. (12.151) . Hence, for large field values 
it is possible to satisfy WMAP normalization and the constraint on r simultaneously for 
much larger values of T„/m, as compared to the quadratic case. 
The non-linearity parameters are given by 



NL 



5Nl_ 
6iV'2' 



9nl 



25 iVf 
54 iV'3 ' 



(3.16) 
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Figure 3: The red (solid) line shows r for the cosine potential and the blue (dotted) line 
for the quadratic potential, respectively, for Ta/m = 10~^ and 10~^. r is seen to decrease 
as T(j/m decreases. 



where 



tot I 



and N': 



tot 



?3 . 

ini 



(3.17) 



dOL *' * del 

Fig. m shows the numerical results of /nl for axion-type (red solid line), quadratic (blue 
dotted line), and quadratic with correction (green dashed line) models. We see that the 
three results agree very well for small 6*^, and begins to differ as 6^ increases. The difference, 
however, is only of order one and hence small. Fig. H] gives a confirmation of the results 
of [12]. The correction case is not really correct at large but we still plot it to show the 
agreement at small 6'*. 

Next we look at the results for qml for the same three models as above. The numerical 
results are plotted in Fig. [5l The red solid line shows the result for the axion-type curvaton, 
and we can see that Qnl is positive and larger than the magnitudes for the other two cases. 
For the quadratic case, we know from the analytic expression of Eq. (12. 9p that it is negative 
and that is shown by the blue dotted line which gives the exact numerical answer. Again 
we know from Eq. ( 12.23^ that for correction case it is positive and that is what we see 
from the green dashed line in the figure. 

From Fig. E] where we have plotted gisn versus fjsn an interesting observation is that 
qnl almost linearly depends on fi\iL for the axion-type curvaton model, namely 

9nl ^ c/jvi, (3.18) 
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Figure 4: Plot of f^i- The red (solid), blue (dotted) and green (dashed) lines correspond 
to axion-type, quadratic and quadratic with correction models, respectively. 
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Figure 5: Plot oi qnl- The red (solid), blue (dotted) and green (dashed) lines correspond 
to axion-type, quadratic and quadratic with a'^ correction models, respectively. 
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Inl 

Figure 6: Plot of g^L vs. f^L for axion-type curvaton. 

where c ~ 732. From Fig. H] and [5l we can see that the formulas for the quadratic with 
correction capture some features of the axion-type curvaton model. Here we try to use 
them to roughly understand why (yfjvL almost linearly depends on /jvl- At the leading 
order, we have 

/i2 ~ /is ~ 1.193(5 = 0.05^2^ (3.19) 

and then 

hL^ 1.25/ fn, gNL^ 0.2161/ fl. (3.20) 
Taking into account the analytic formula for in Eq. fl2.19p . we obtain 

fNL ^ w-l-J'-e;', (3.21) 

M'^ T 

g^L ^ l3.A-l—e;\ (3.22) 
Therefore, the coefficient c is roughly given by 




(3.23) 



Using f /Mp = 5 X lO"'^ and T„/m = 10"^, c ~ 536 which is roughly the same as that for 
the axion-type curvaton model. A general discussion for the curvaton model with a small 
correction to the quadratic potential is given in the Appendix. 
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3.2 Mixed scenario 



In the previous subsection, where we assumed that the total primordial power spectrum 
is generated by the curvaton only, we saw that, having fixed f /Mp = 5 x 10~^, the value 
of T(j/m = 10"^ is too large for r to satisfy observational constraints. We know from 
Eq. fl2.15p that r scales linearly with T„/m for the quadratic case. From Fig. [3] it seems 
that for the axion-type case also r scales linearly with T„/m, though it is not obvious 
from the equations. Hence, this value was enough to demonstrate the physics. 

In this subsection we consider the scenario where the final curvature perturbations are 
contributed by both infiaton and curvaton fiuctuations. This scenario has been invoked 
in [8] to naturally achieve a red-tilted primordial power spectrum in the curvaton model. 
It is convenient to introduce a new parameter /3, 

/3 = n./n. (3-24) 

where is the total primordial power spectrum including the contribution from infiaton, 
and hence (3 G [0, 1]. Then the effective non-Gaussianity parameters becomes 

rJi=(3'fNL, g^i = P'9NL. (3.25) 

See [8,11] for detail. The same result of f^l in the mixed scenario is also obtained in [13]. 
If 9nl = c/nl, 

9%^l = PcfJi (3.26) 
where g^^[ is suppressed by a small factor j3 with respect to f^l- In [11], the author 
found 

gNL = cfliL (3-27) 

for the curvaton model in which the curvaton has a polynomial potential and the quadratic 
term is subdominant during inflation. In this case, 

= -p{rji)\ (3.28) 

which is enhanced by a factor l/[3. A large qnl is expected in the curvaton model. In 
general, > ^{fjJlY multi-fleld case [19]. In particular, in the mixed scenario, 

T^l is enhanced by a factor 1//3, namely 

r'Jl = §-^irJlr- (3.29) 
For the axion-type curvaton model /? can be expressed as 

- ^ = [K'j^. (3.30) 
14 
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Figure 7: Plot of /^/. The red (solid), blue (dotted) and green (dashed) lines correspond 
to axion-type, quadratic and quadratic with correction models, respectively. 

If r < 8f'^/{N^Mp), then jS < 1, which implies that the quantum fluctuation of the 
curvaton only contributes a part of the total primordial curvature perturbation. In Figs. U\ 
Eland [HI we have plotted /^^, g'jj[ and r^^/, respectively. The parameter values we have 
used for these plots are f/Mp = 5 x 10^'^, T^j/m = 10"^ and r = 0.1. We see that the 
values are suppressed by factor of for fiy^ and by P'^ for g^L and tj^l. 

4 Conclusion 

The curvaton model is one of the few models for the generation of primordial density 
perturbations which predict large non-Gaussianity and is potentially testable in the near 
future in upcoming experiments such as Planck. Specific models that have been well 
studied and for which analytic approximate results can be obtained are the simplest ones 
such as a polynomial potential [11]. In this paper we solved the evolution equations for 
the quadratic potential to see whether the analytic approximations are really accurate. 
We found that the analytic expression underestimates the exact answer by about factor 
of two. The reason is that the axion does not exactly behave as dust-like matter when 
it oscillates about the minimum of its potential. The friction term in the axion equation 
of motion, though small, is not negligible and makes the energy density decrease with 
respect to the scale factor as a~^~^'^, where e is small and positive. 
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Figure 8: Plot oi g]^[. The red (solid), blue (dotted) and green (dashed) lines correspond 
to axion-type, quadratic and quadratic with correction models, respectively. 




Figure 9: Plot of t'^J[- The red (solid), blue (dotted) and green (dashed) lines correspond 
to axion-type, quadratic and quadratic with correction models, respectively. 
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Our main motivation in this paper was to understand the non-hnearity parameter 
qnl in the axion-type curvaton model. To achieve this goal we used the SN formalism to 
solve for the two non-linearity parameters. The axion type model is very hard to solve 
analytically for large field values. In addition, computing (^tvl numerically is much harder 
than /tvl since it contains a third order derivative as against the second order derivative 
of I'nl- We have evaluated the tensor-scalar ratio r using WMAP normalization and 
used it to compute Jnl and Qnl- We compare them with results for the curvaton with 
quadratic potential. We found that r is much smaller for the axion-type case compared to 
the quadratic potential when a* ~ /. The difference between them becomes increasingly 
more pronounced at larger field values, qnl is found to be positive, unlike the quadratic 
case where it is negative, and the amplitude of g^L is much larger. We also found that 
there is a nearly linear scaling between gjsfL and Jnl, with small deviation from linearity 
at large field values. Lastly, we considered a mixed scenario where both the inflaton and 
the curvaton contribute to the primordial power spectrum. We have found that g^i and 
Jnl are greatly suppressed with respect to the case where only the curvaton contributes. 
Moreover, g'^J[ is suppressed with respect to /^^, whereas, r'^l is enhanced and hence it 
would be inconsistent to ignore the contribution to the overall primordial non-Gaussianity 
from t'^I- 

We conclude with a few remarks on the trispectrum. Usually the trispectrum is 
characterized by two parameters tnl and gNL- While J'nl is well studied for a variety of 
curvaton models in the literature, the trispectrum is not as well studied. One reason is 
the technical difficulty of having to deal with third order perturbations. Several recent 
works have shown that tnl and g^L can be as important as, if not more than, f^i- 
On the observation side much work has been devoted to constraining the value of Jnl 
using the CMB observations. It makes sense to ignore t^l and gNL their contribution 
is small. However, we need to take them into account if the trispectrum is large since 
they are independent parameters and have the potential to discriminate between different 
curvaton models. It would be very interesting to study the implications of large tnl and 
g^L on the CMBR and we are working along those lines. 
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A qnl vs. /iVL in the curvaton model with small de- 
viation from quadratic potential 



The non-Gaussianity for the curvaton model with small deviation from quadratic potential 
is investigated in [10]. In this appendix we extend the previous discussion. Assume the 
curvaton potential is given by 

y(a) = imV + Am^f-)", (A.l) 

where n > 2. U n < 2, the correction term becomes dominant around a = 0. The 
correction in the equation of motion is small if |s| <^ 1/n, where 



s = Xr- . (A.2) 
Vm/ 

The general result is given in [10]. Here A = — m^/(24/^) and s = —6 for the axion-type 
curvaton model. For case of the small correction, Jnl and qnl take the form 

fNL ^ 7^, (A.3) 



25 , , s 25 



9nl ^ -777^(^^)72- - (A-4) 



where 

v{n) = -n^{n - \){n - 2)2'^r(5/4)'"-V 

^i(-) /' Jr'(x)Fi(x)x^rfx-Fi(-) J1{x)x^dx , (A.5) 
■i 2 Jq 4 4 t 2 i 

which shows up in Fig. [TOl From Eq. flA.4l) . gj^n becomes positive if 

s < -^Jd. (A.6) 

v{n) 

If Qnl is dominated by the second term in Eq. (1A.4I) . g^L — —^Inl- Now we focus on 
the case of |s| ^ Af^i/viji), namely 

25 s 

9nl ^ -^^H jT- (^•'^) 

Using Eq. fl2.14l) . we have 



200 .,Kr. /a,x-6 

9nl ^ — —v{n)X^— — . (A.9) 
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Figure 10: Plot of v{n) as a function of n. 



Therefore, 

^^,.^c/J7/^ (A.io) 

where 

For n = 2, c = v{n) = and qnl still linearly depends on f^i. In the mixed scenario, 

9''Jl = f3''-'c{rJir''^'. (A.12) 

The second order non-Gaussianity parameter is enhanced in the mixed scenario if n < 3. 
However, for 2 < n < 3, v{n) < 1, and hence a large absolute value of g^L is quite 
unlikely. 
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